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Abstract 

We adduce the results of the condensate fraction calculation for 
liquid helium-4. The method is derived from the first principles and 
involves minimum assumptions. The only experimental quantity we 
need in our calculations is the static structure factor which is easily 
measurable. We use the approximation in which expressions contain 
one summation in the wave vector space (or one integration it the ra- 
dius vector space) in order to demonstrate the validity of our method. 
The computed values of the condensate fraction lie within 8.8 -r- 14% 
depending on the model potential for the short-range interaction. This 
result agrees with recent experimental measurements and numerical 
estimations. 
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1 Introduction 



The phenomenon of Bose-Einstein condensation and superfluidity in Bose- 
systems is a subject of great interest. Superfluidity in liquid helium-4 can 
be easily observed now but the measurements on the condensate fraction 
remain a difficult problem because of the strong interactions in this fluid. 
Researchers, both experimentalists and theorists, now generally concur that 
~ 10% of atoms are in the lowest (zero momentum) state at zero temperature. 
The relative number of atoms with zero momentum is called condensate 
fraction, here we denote it as /. 

The problem of the calculation of condensate fraction in liquid helium-4 
was first considered in the classical paper by Penrose and Onsager [3|. The 
authors used a very rough approximation of hard spheres for the ground-state 
wavefunction and calculated this quantity as approximately 0.08 or 8%. 

During recent twenty years different methods for measurements and theo- 
retical calculations of the condensate fractions were used. Neutron scattering 
is used to draw the information on the momentum distribution and hence 
on the condensate fraction. Sears et al |2j obtained the value 13.9%. Later, 
similar value of 13.3% was obtained by Sears |3j using the data on the tem- 
perature variations of the average kinetic energy. 

Numerical study based on the variational method were carried out in the 
series of works by Manousakis, Pandharipande et al 0, j|, § by means of 
Jastrow wavefunction with three-body corrections. Their results vary within 
8.2 -=-10.3%. 

The numerical estimations of the condensate fractions were made also in 
the series of works by Vakarchuk et al and Vakarchuk |7|, || [J leading to the 
results of 3.7 -j- 8.8%. The methods used in these papers are based on the 
mean-spherical approximation for the structure functions of liquid helium-4 
at zero temperature direct quantum-mechanical computations from the 
first principles [§, |J . 

The value of approximately 10% was obtained for / by Sokol and col- 
laborators by means of deep-inelastic neutron scattering at high-momentum 
transfer §i§. 

Monte Carlo (MC) methods were also used to study the question of the 



condensate fraction. The result of Whitlock and Panoff [13] is ~ 9% (Green 



Function MC). Ceperley and Pollock [|14]] obtained value ~ 7% using Path 
Integral MC technics at 1.18 K. Also, Moroni et al [15 gives 7.2% from 
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Diffusion MC simulations. 

Approximately in the same time the experimental result of ~ 6% was 
obtained by Azuah et al |jrB| . 

Series of papers by Mayers et al Jl7|, [18[ was devoted to both calculation of 
the condensate fraction based on phenomenological assumptions [|17]] (result 



is 9.9%) and measurements on the high-energy scattering [|18j (/ = 15 ± 4% 
at 1.3 K). 

Recent high precision measurements of the dynamic structure factor al- 



lowed Glyde et al [|19] to derive the information on the condensate frac- 
tion at zero temperature and its temperature dependence. They report 
/ = 7.25 ± 0.75% at K. 

A semi-phenomenological method for extracting the condensate fraction 



data was proposed recently by Rinat and Taragin [20|. Authors obtained / at 



different temperatures and extrapolated their results to K as / = 9.0±0.3%. 

The aim of this work is to show a possibility of calculation of the conden- 
sate fraction without drawing any special huge computational efforts. 

The method we use was worked out in ||. The condensate fraction is 
calculated directly from the single-particle distribution function Fi(R) as its 
long-range limit: 

f = N /N = lim F 1 (R) 

R^oo 

where Nq is the number of particles in zero momentum state and N is the 
total number of particles in the system. 

Our technics do not involve any phenomenological assumptions. For the 
computations we use as input only one experimental quantity, the static 
structure factor of helium at zero temperature. It is the advantage comparing 
to those methods for which such quantities as the dynamic structure factor 
is necessary. Its precise measurement is much more complicated problem 
comparing with that of the static structure factor. 

The method is fully controlled, i. e., we always know which effects are 
neglected, and it is not difficult to take them into account in the next ap- 
proximation (in principle). Here we present only the first approximation 
(RPA) but the expressions are easily extendable in order to take subtler 
effects like many-particle interactions. 

The only key problem arising while extending our method further is the 
knowledge of the static three-particle structure factor. Unlike usual structure 
factor, it is not an easily measurable quantity already. One of the ways for 
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its calculation we see in such a scheme: expressions for the three-particle dis- 
tribution function F3 (beyond superposition approximation) containing pair 
distribution function — > 6-dimensional Fourier transformation of F 3 leading 
to the three-particle structure factor. A separate paper will be devoted to 
this problem. 

We consider three different approaches for the short-range repulsive in- 
teractions. In the first one, they are not taken it into account explicitly at 
all. In the second approach this repulsion is modeled by Meyer's function 
e (A/R) n — \ r [ s the radius- vector absolute value. We chose A = 2.1 A and 
n = 12. This potential is referred here as 'almost hard spheres' (AHS) [pl[ |. 
The third model is the potential of hard spheres (HS) with diameter 2.1 A. 



2 Calculating procedure 

Principal calculating procedure consists of the following steps: 

1. Input: Experimental data on the structure factor [^] (should be con- 
verted into the zero temperature [23]); 



2. Input: Model for the short-range interactions (HS, AHS or none) [21 

3. Output: Structure factor of the model system corresponding to the 
short-range repulsion; 

4. Output: Fourier transformation of the interatomic potential (effective) 

5. Output: Pair distribution function F 2 (R) f22|| ; 

6. Result: Condensate fraction. 

In this work, we propose a simple method for the calculation of the con- 
densed fraction using formulae obtained earlier in ||: 

/ = exp(/ 1A + A 1 J + ...), (1) 
IlA = "Ia?^^ -1 ) 2 /^' 

Ai J = gjdR [2h*(R) - h(R) + h 2 (R)/A 



h(R) = F 2 (R)-l, h*(R) = JF 2 (R) - 1 
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dots denote terms having more than one integration over wave vector or 
coordinate, F 2 (R) is the pair distribution function (PDF), N is the number 
of atoms, and 

q is helium density, g = 0.02185 A~ 3 , m is the mass of helium atom, m = 
4.0026 a. m. u., q is the wave vector. 

The expression for PDF within the accepted approach was obtained in [HI: 



I N ^ 1 



^)=e X p^E^ J e- + ...J, (3) 




while the self-consistent description of short- and long-range correlations in 
liquid helium-4 results in such a result for the function a 2 pif : 

(4) 

with Sg tm being the structure factor of the model system (it equals 1 if the 
short-range repulsion is not given explicitly), superscript 's. r.' corresponds 
to 'short-range'. S g is the experimental structure factor (at K). These 
expressions contain v q being the Fourier image of the interatomic potential 
in helium that has an effective nature. It is obtained from the first princi- 
ples using the collective variables formalism in the Schrodinger equation, as 
described in p5|. 



On the step (|4|) of the calculation procedure we encountered problems 
while considering the HS short-range potential. The reason is clear: a weak 
damping of the structure functions at large values of wave vector or ra- 
dius vector requires much more careful numerical computations comparing 
to those in two other cases. 

Some usual problems appear on the step (§) of the calculating procedure 
where we smoothed F 2 {R) calculated from (|3|) and truncated it at the dis- 
tances R < 2.2 A. We put F 2 (R < 2.2 A) = in order to obtain a correct 
PDF at small distances |25|| . We present the calculated PDFs in comparison 
with the experimental (being more precise, derived from the measurements 
on the liquid structure factor) one [^| in the Fig. [TJ. 
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Figure 1: Pair distribution function. The dashed line corresponds to the 
experimental data by Robkoff and Hallock [25] at 1.38 K. The dotted line 
represents the smoothed data in the case when no explicit short-range repul- 
sion is given. The dashed-dotted and solid lines are PDFs for HS and AHS 
short-range potentials respectively. 



3 Results and discussion 

We have calculated the condensate fraction at zero temperature using differ- 
ent models for the short-range repulsive part of the potential. 
The calculated values of / are: 

• No explicit short-range potential / = 14%; 

• AHS short-range potential / = 11%; 

• HS short-range potential / = 8.8%. 



We expect the real value being about 10 30% less than the calculated 
one due to the contribution of the higher-order terms, as described in M. 
But in the first approximation we consider it to be a good accuracy for / in 
liquid helium-4. 

Our result for the condensate fraction is within the range of theoretical, 
numerical and experimental estimations mentioned here in the introductory 
part of the paper. Since our expressions are derived from the first princi- 
ples, one can use the number for / calculated in this way as a test for the 
condensate fraction information indirectly extracted from the experimental 
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measurements. This becomes possible after the next one and two approxi- 
mations are computed. 

We found that the short-range repulsive part of the potential has a sig- 
nificant influence on the results for the condensate fraction. Namely, when 
the short-range interactions are included into the consideration explicitly, the 
results appear to be essentially dependent on the "hardness" of the core. A 
proper model for this interaction will be found when more terms are taken 
into account in the expressions (|l|), (J3|). The reason is that in principle the 
summation of the whole series corresponds to the correct calculation of the 
short-range interactions pi. On the other hand, the condensate fraction 



in liquid helium-4 is a very sensitive (with respect to the order of approx- 
imation) quantity. Thus, we can expect that having found such a model 
we obtain a possibility to use simple RPA-like expressions instead of higher- 
order ones. By now, the AHS (n = 12, A = 2.1 A) potential seems to be the 
most suitable for this purpose. 
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